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The semileptonic process, B — > n lu, is studied via full QCD lattice simulations. We use un- 
quenched gauge configurations generated by the MILC collaboration. These include the effect of 
vacuum polarization from three quark flavors: the s quark and two very light flavors (u/d) of vari- 
able mass allowing extrapolations to the physical chiral limit. We employ Nonrelativistic QCD 
to simulate the b quark and a highly improved staggered quark action for the light sea and va- 
lence quarks. We calculate the form factors f+(q 2 ) and fo(q 2 ) in the chiral limit for the range 

f6GeV 2 < q 2 < and obtain f*™?* v2 [ dT/dq 2 ] dq 2 / \V ub \ 2 = 1.46(35) ps" 1 . Combining this 

with a preliminary average by the Heavy Flavor Averaging Group (HFAG'05) of recent branching 
fraction data for exclusive B semileptonic decays from the BaBar, Belle and CLEO collaborations, 
leads to \V ub \ = 4.22(30)(51) x 10 -3 . 

PLEASE NOTE APPENDIX B with an ERRATUM, to appear in Physical Review D, to the pub- 
lished version of this e-print (Phys.Rev.D 73, 074502 (2006)). Results for the form factor f+(q 2 ) in 
the chiral limit have changed significantly. The last two sentences in this abstract should now read; 
"We calculate the form factors f+(q 2 ) and /o(<? 2 ) in the chiral limit for the range 16GeV 2 < 

2 

q 2 < q^ax and obtain § x ^ e X y 2 [ dT/dq 2 ] dq 2 / \V u b\ 2 = 2.07(57) ps _1 . Combining this with a pre- 
liminary average by the Heavy Flavor Averaging Group (HFAG'05) of recent branching fraction 
data for exclusive B semileptonic decays from the BaBar, Belle and CLEO collaborations, leads to 
\V ub \ = 3.55(25)(50) x 10" 3 ." 

PACS numbers: 12.38.Gc, 13.20.He 



I. INTRODUCTION 

A major achievement of the B factories in recent years 
has been the observation of CP violation in the neutral 
B system [l|, The emphasis since then has been on 
overconstraining the unitarity triangle and checking for 
consistency with, or deviations from, the three family 
Standard Model. The goal is to independently measure 
not only the three angles but also the lengths of the sides 
of the triangle and thereby determine (p, rf) , the apex of 
the unitarity triangle, in as many ways as possible. In 
order to fix the sides of the unitarity triangle, the mag- 
nitudes of several Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements are required and the accuracy is cur- 
rently limited mainly by theoretical uncertainties in two 
of them, \Vtd\ and \V u b\. Theory input for these quan- 
tities involves hadronic matrix elements of several elec- 
troweak operators and these in turn require good control 
over nonperturbative QCD. This article reports on sig- 
nificant recent progress in lattice QCD determinations of 
form factors relevant for the CKM matrix element \V u b\. 

\V u b\ can be determined from studies of either inclusive 
or exclusive B meson semileptonic decays. The first de- 
terminations relied on inclusive measurements. However, 
recent impressive progress in measurements of branching 
fractions for exclusive decays by CLEO Q, Belle [1, [ff 
and BaBar 0, 0, [1] have started to make exclusive deter- 
minations competitive [tj U3] ■ In either approach errors 



are now dominated by theory errors and which method 
will eventually win out depends on how well the theo- 
retical uncertainties in shape functions (for the inclusive 
approach) or form factors (in the case of the exclusive 
approach) can be brought under control. Lattice QCD 
provides a first principles nonperturbative QCD method 
for calculating form factors in exclusive semileptonic de- 
cays. The first lattice calculations were carried out in the 
quenched approximation that ignored vacuum polariza- 
tion effects [lj [H [H H For some period these 
pioneering results were the only ones available and ex- 
perimentalists, e.g. the CLEO collaboration, used them 
to extract some of the earliest exclusive \V u b\ results [3[. 
In the summer of 2004 the first preliminary unquenched 
results for the B — ► irlv form factors we re p resented 
by the Fermilab/MILC [16] and HPQCD [If collabo- 
rations. These calculations employed the MILC collabo- 
ration Nf = 2+1 unquenched configurations, the most 
realistic gauge configurations to date with vacuum po- 
larization from 2 flavors of very light quarks and from 
strange quarks [18j . Furthermore, the good chiral prop- 
erties of the improved staggered quark action used for the 
light sea and valence quarks, allowed for investigations 
much closer to the chiral limit than in the earlier calcula- 
tions. Hence, references [l?} constitute a major step 
forward in lattice QCD determinations of semileptonic 
form factors. These preliminary results have been incor- 
porated by Belle and BaBar into their recent B semilep- 
tonic analysis 0, H, IE & HI and used by other theorists 
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in their extractions of \V u b\ fl9l I20I l2ll| . 

References [IH E3] are part of a growing list of recent 
lattice calculations that use the MILC unquenched con- 
figurations. The creation of these configurations became 
feasible on present day computers due to the develop- 
ment of highly improved staggered light quark actions 
[22I ] . There is one well known drawback of staggered ac- 
tions, namely that each flavor comes in four types, called 
"tastes" . To simulate just one taste of sea quark per fla- 
vor a fourth root of the quark determinant is used and 
the validity of this procedure has not yet been rigorously 
proven. All tests undertaken to date to study the fourth 
root procedure, however, have led to encouraging results 
[23I ]. There are no indications of problems or any de- 
viations from continuum QCD behavior beyond small 
and expected discretisation errors that can be system- 
atically improved upon. A recent review of the fourth 
root issue is given in (24j|. The MILC unquenched con- 
figurations and the heavy and light quark actions cur- 
rently in use have also been tested repeatedly by cal- 
culating a large set of well measured quantities such as 
light hadron spectroscopy, light meson decay constants, 
quarkonium and B meson spectroscopy. Agreement be- 
tween experiment and lattice calculations is found to 
be excellent within the few percent errors in the lat- 
tice results [H, H^, H3, HI, HI] - More impressively, these 
same gauge configurations have led to some recent pre- 
dictions from Lattice QCD that have been confirmed sub- 
sequently by experiment [2(3, EHl H3] • Most relevant for B 
semileptonic decays, is the successful calculation by the 
Fermilab/MILC collaboration of D meson semileptonic 
form factors and the agreement of their q 2 dependence 
with experiment [13]. All this should boost confidence 
in results for B meson semileptonic decays that are now 
emerging from Lattice QCD calculations. 

In this article we make several further improvements 
on the preliminary calculations of ref. [13] and finalize our 
form factor results. Some of these improvements have 
been reported on in [33| . Among the improvements, we 
have now included all the dimension 4 current corrections 
through 0{a s /M) and 0{aa s ) to the temporal and spa- 
tial components of the heavy-light electroweak currents 
relevant for B — > n Iv semileptonic decays. Our previ- 
ous results included currents only at lowest order in 1 /M 
and through 0(a s ). We now have simulation data from 
several MILC ensembles, whereas in ref.fl7} results were 
at a single fixed light sea quark mass and only the va- 
len ce q uark mass was varied. Another development since 
ref. [13] is that we now employ the formulas of staggered 
chiral perturbation theory [34 , l35l |36( simultaneously to 
both partially quenched and full QCD results to extrap- 
olate in the light quark mass to the physical pion. Previ- 
ously we used linear chiral extrapolations to our partially 
quenched data. 

So, the list of improvements in our semileptonic de- 
cay studies since the work of ref. [13] is quite substantial. 
Nevertheless, we will see that changes in the final re- 
sults for the form factors f+(q 2 ) and fo(q 2 ) are almost 



negligible. The individual contribution from each higher 
dimension current correction is found to be small, mainly 
due to small one-loop matching coefficients. There is also 
some cancellation between the many terms. The differ- 
ence between linear chiral extrapolations of our previous 
partially quenched data and results from staggered chiral 
perturbation theory fits to the new full data set also turns 
out to be a small effect. With all the improvements now 
in place, our semileptonic form factor calculations are at 
the same level as recent B meson decay constant deter- 
minations by the HPQCD collaboration [53, - Those 
latter calculations are crucial for determinations of the 
CKM matrix elements \V t d\ and |V ts |. 

In the next section we provide some details of the sim- 
ulation parameters and of the lattice actions employed. 
We also summarize formulas for the relevant form factors. 
In section III. we discuss matching of the lattice heavy- 
light currents used in our simulations to their continuum 
QCD counter parts. One-loop matching coefficients for 
the temporal components Vo and A were published in 
ref. [39( and have been used already in our decay constant 
determinations [33, HH ■ Matching of the spatial compo- 
nents Vk was completed as part of the current project. 
Section IV. describes how we extract the form factors /|| 
and f± from numerical simulations. The more commonly 
used form factors /+ and fo can be expressed as simple 
linear combinations of /|| and fx- Section V. focuses on 
chiral extrapolations of f» and fx- In section VI. we 
summarize our final results for f+{q 2 ) and fo{q 2 ) in the 
physical chiral limit. We also present Tables of partially 
integrated differential decay rates (divided by |Kb| 2 ) in 
several q 2 bins. In section VII. we combine the results of 
section VI. with experimental data on B — > tt Iv branch- 
ing fractions to estimate \V U b\- We then conclude with a 
summary section. 



II. SIMULATION DETAILS AND FORM 
FACTOR FORMULAS 

Most of our simulations were performed on the MILC 
collaboration "coarse" ensembles with lattice spacing 
around 0.1 2fm and with light sea quark masses in units of 
the strange quark mass, mj/m s , ranging between 0.125 
and 0.5. We have also carried out some checks using 
MILC lattices with finer lattice spacings. We refer to the 
original papers by the MILC collaboration for details of 
how their lattices were generated [l8|. Some specifica- 
tions for the ensembles are given in Table I. The highly 
improved staggered action, the AsqTad action, is used 
for both the sea and valence light quarks. The lead- 
ing discretisation errors in this formalism are 0(a 2 a s ). 
To simulate the b quark inside the B meson we use the 
same highly improved nonrelativistic (NRQCD) action 
employed in recent studies of the T system on the same 
MILC ensembles 0, [4(J ■ In terms of the two-component 
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TABLE I: Simulation Details, m/ (ra q ) denotes sea (valence) 
quark masses and m s is the strange quark mass. A f means 
the lattice spacing was determined through n. For all other 
ensembles the T 2S-1S splitting was used, uo = [plaq] 1 ^ 4 is 
the link variable that enters into the MILC collaboration's 
convention for normalisation of quark masses. 
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where n is a stability parameter introduced to control 
high momentum modes of the b propagator [28|. Hq is 
the nonrelativistic kinetic energy operator, 
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and SH includes relativistic and finitc-lattice-spacing cor- 
rections, 



C4 on B semileptonic form factor calculations to be less 
than 1%. 

Just as in continuum QCD, our lattice actions include 
a small number of parameters that can only be fixed via 
experimental input. These are the bare quark mass pa- 
rameters and the scale (or coupling). For the present case 
of heavy-light simulations the action parameters have al- 
ready been fixed for us via simulations of quarkonium 
and light quark systems using the same MILC configu- 
rations. We use lattice spacings determined by the T 
2S-1S splitting [Hj]. For two of the MILC ensembles (de- 
noted by a f in Table I), T simulations have not yet been 
carried out. On the other hand, the MILC collaboration 
has measured the heavy quark potenti al p arameter r\ja 
(in lattice units) for these ensembles [26|. In reference 
[28], using MILC ensembles on which both r\ja and the 
T 2S-1S splittings were calculated, we could use the ex- 
perimental 2S-1S splitting to determine r\ — 0.321 (5)fm. 
This physical value for r\ was combined with the MILC 
collaboration r\ /a to fix the scale for the above two en- 
sembles. The b quark mass is also fixed by our T studies. 
Studies of pions and kaons have fixed the u and d masses 
(which we take to be equal to each other) and the s quark 
mass respectively [III [27j ■ Hence by the time one gets to 
the B system there are no adjustable parameters left in 
our QCD action. 

To study the process B — > irlv, one needs to evaluate 
the matrix clement of the charged electroweak current 
between the B and the ir states, (7r|(V — A)^\B). Only 
the vector current contributes to the pseudoscalar-to- 
pseudoscalar amplitude and the matrix element can be 
written in terms of two form factors f+(q 2 ) and fo(q 2 )- 
These depend only on the square of the momentum trans- 
fered between the B and the 
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If one neglects the mass of the charged lepton in the 
final state, only the form factor f+(q 2 ) contributes to 
the decay rate T(B — * nil/). Nevertheless, it is useful to 
keep track of the form factor fo(q 2 ) as well since, as we 
shall see, it helps in our interpolations and extrapolations 
of simulation data. In our data analysis another pair of 
form factors, f» and f±, turn out to be more convenient. 



Mq is the bare 6-quark mass, A' 2 ' the lattice Laplacian, 
V the symmetric lattice derivative and A^ 4 ) the lattice 
discretization of the continuum J~). Df. Expressions for 
the improved E and B fields are given in [28| . All deriva- 
tives are tadpole improved. As in [28| we set the c,-'s to 
their tree level value, c; = 1. These coefficients will be 
modified by radiative corrections at higher order. For 
heavy-light systems the most important radiative correc- 
tions will be to C4. All other Cj's are multiplied by addi- 
tional factors of l/M or a. Based on the discussions in 
12811 we estimate the effect from radiative corrections to 



(n{p.)\V»\B{p B )) = ^2M B K/|| + p^/jj, (5) 



with 



u 

Pb 
Mb 



rt-{p.-v)v». (6) 



In the B rest frame (in this article we only consider B 
mesons decaying at rest) the temporal and spatial parts 
of ([5]) become, 



(7r\V°\B) = v^Mfl/n 
(n\V k \B) = V^pt/i 



(7) 
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Hence, one sees that one can separately determine /|| or 
f± simply by looking at either the temporal or spatial 
component of V^. These two form factors have the ad- 
ditional advantage that they have simpler HQET scaling 
properties and chiral perturbation theory is carried out in 
terms of them rather than for /+ and Jo- After carrying 
out the chiral extrapolations for f» and f±, we convert 
back to obtain /+ and /o for the physical theory using, 
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where E n is the pion energy in the B rest frame. From 
these formulas one sees that / + will be dominated by f±, 
i.e. by the matrix element of V k , and fo by /|| or the 
matrix element of Vq. 

Our goal is to evaluate the hadronic matrix elements 
(7r|U°|i?) and (7r|U fc |-B) via lattice simulations. There are 
several steps in the calculation. First, one must relate the 
continuum electroweak currents, V° and V k , to lattice 
operators written in terms of the heavy and light quark 
fields in our lattice actions. In the second step the matrix 
elements of these lattice current operators must be eval- 
uated numerically and the relevant amplitude, i.e. the 
matrix element between the ground state B meson and 
the ground state pion with appropriate momenta must 
be extracted. This will give us, via eqns.([7]), the form 
factors /|| and f± as functions of the light quark mass 
and the pion momentum. Finally in step 3 these numeri- 
cal results must be extrapolated to the physical pion. In 
the next three sections we describe each of these three 
steps in turn. 



III. MATCHING OF HEAVY-LIGHT 
CURRENTS 

Matching of heavy-light currents between continuum 
QCD and a lattice effective theory with two-component 
nonrelativistic heavy quark fields $ and four-component 
light quarks q(x) is discussed in ref. [4l[. Since stag- 
gered light quarks can be written in terms of four- 
component "naive" AsqTad quark fields the formalism 
developed there carries over unchanged to the present 
calculation. Introducing also a four-component notation 
for the heavy field, Q(x) = ($, 0), one finds that through 
0(a s Aqcd/M , a s /(aM) , a s aA QC D) the following 
current operators in the effective theory are required. 

temporal : 



4°\-) 



q(x) T Q(x), 
^(x)T Ql .VQ(xl 

l T q{x)-f^ l0 r Q(x). 



2M 



(10) 



spatial : 



4 0) (x) = q(x)T k Q(x), 
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q{x) V fe Q(x), 



J, 
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0*0 = ^n~9(a;) v" fe Q(x), 



2M, 
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where T M can be either 7^ or 757^, and Mq is the bare 
heavy quark mass in the NRQCD action. One sees that 
there are two dimension 4 current corrections to the tem- 
poral component and four such corrections to the spatial 
components. To the order that we are working, one has 



(V ) = (l + a s p^)(4 0) ) + 

(l + asP ^)(J^ sub )+a s p^(4 2) ) (12) 
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We introduce the combination J, 



(l),sub 



J, 



(i) 



u s Cio,/i^m • This subtracts out power law contributions 
to the matrix elements of the higher dimension operator 
4 1] through 0(a s /{aM)) 0. enters the matching 
already at tree level and after the subtraction one is left 
with the physical 0(Aqcd /M) contribution that is a rel- 
ativistic correction to the leading order term. Power law 
subtractions of the other dimension 4 current corrections 
come in as 0(a 2 s /(aM)) effects and are only partially in- 

(i) 

eluded here. The one-loop coefficients p)l and £10,^ for 
(j, = are given in ref.[39|]. The results for p = k have 
not been published before and are summarized in Table 
II. In ref.[17| only the contributions from the first lines 
in eqns. (fl"2] ) and (fl~3|) were taken into account, i.e. Jq ' 
and matched through C(a s ). 

As mentioned in the introduction, the effects of all the 
dimension 4 current corrections turn out to be very small. 

In Fig.l we show results for (4 )/(4 ) l° r one °^ our 
ensembles with and without the power law subtraction. 

One sees that although the unsubtracted ^o°^) is 

at the - 6% level, the physical ( J^' 15 " 6 )/^) is < 1%- 
In Figs. 2 - 4 we give further examples of (J k ^) / (J^) 
for j > 1. These get multiplied by factors of (p k ^a s ) in 
eq.©. Using a s w a v (2/a) = 0.32 [H and Table II, 
one finds pa s factors between 0.01 and 0.11, which leads 
to contributions from higher order currents that are at 
most 1%. For instance the largest current correction is 
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0.063(3) 


0.073 



TABLE II: Matching coefficients for the spatial currents 14. 
Where errors are not indicated explicitly, they are of order one 
or less in the last digit. aMo is the bare heavy quark mass in 
lattice units and n a parameter in the NRQCD action. The 
three selected values for aMo correspond to the b quark on 
the MILC extra-coarse, coarse and fine lattices respectively 
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FIG. 1: The ratio (J (1) ) / (J (0) ) for one ensemble (uoamf — 
uoam q = 0.01) versus the pion energy E n . The lower points 
are before power law subtraction and the upper points after 



power law subtraction (i.e. (Jq ' SIii ') / (Jq U> ) ) 
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FIG. 3: Same as Fig.2 for (jf V^jT) 
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FIG. 4: Same as Fig.2 for {jf*) / {J^) 
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FIG. 2: The ratio (4 ; )/<4 ) for two ensembles versus the 
pion energy E-k. Squares are for uoamf — 0.01 and circles for 
uoamf = 0.02. 



4 4) (Fig.4), but p£> a s = -0.029 a s = -0.0009, and the 
contribution from this current is negligible. In comparing 
Figs. 1-4 one sees that the size of matrix elements grows 



with the pion energy for J^. 2 ' and J^' and seems 
much less sensitive to E v for the other two currents. This 



reflects the fact that jj?^ and J^ 1 ' have derivatives acting 
on the light quark field that is part of the final state pion 
and therefore knows about its momentum. 



(4) 



IV. SIMULATION RESULTS FOR FORM 
FACTORS / AND f ± 

The starting point for calculations of the hadronic ma- 
trix elements (n\J^\B) is the 3-point correlator, 



C^(p^p B ,t,T B ) = 
££>.(0)J«(z,i)4(£-T B ))e 



iPB-V e i(Pw-PB)-z 



(14) 
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where <I>b and $,r are interpolating operators for the B 
meson and the pion respectively. All results here have 
the B meson three momentum, ps, set equal to zero. 
For simplicity, the pion operator (Ev was always placed 
at the origin. The B meson was then created at time 

(k) 

slice —Tb and the electroweak current, , that con- 
verts the b quark into an u quark was inserted at times 
> t > —Tb- We have also simulated the time- reversed 
process, which then has the electroweak current inserted 
between +Tb >t>0 and <f>s acting on time slice +Tb- 
By looking at both forward and time-reversed processes 
and verifying that they lead to consistent results (within 
statistical errors), we were able to increase statistics and 
at the same time provide some check on our codes. For 
most of our simulations we used Tb/cl — 16 on the coarse 
MILC lattices and Ts/a — 24 on the fine lattices. On 
one of the coarse lattices we also ran with Tb /a = 20 and 
verified that results for form factors were independent of 
Tb- Making Tb too large is not helpful since statistical 
errors grow with Tb ■ On the other hand making Tb too 
small limits the number of data points available and gives 
us less flexibility in our fits. 

In constructing the interpolating operators and ^ 
we have found it convenient, just as in the currents of 
cqs.(|f Oj) and pip , to work with four component naive 
fields q(x) rather than one component staggered fields 
X(x). Hence in eq. p4[) we use, 
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The relation between naive and staggered fermion prop- 
agators is given by [44] , 

G q (x,y)=n(x)rt(y)G x (x,y) (16) 



with 

3 

fj,=0 

In our simulations we first calculate staggered propaga- 
tors G Xl since they are cheaper, and then convert to naive 
propagators G q using eq- fTB]) before evaluating 3-point 
correlators. The naive AsqTad theory has 16 tastes of 
quark per flavor and hence one could form 16 different 
heavy-light pseudoscalar boundstates. However, as dis- 
cussed in ref. [H| these exactly degenerate 16 B mesons do 

not mix and the 2-point correlator {&b{%)®b(v)) for in- 
stance, receives contributions from only one of these pos- 
sible B mesons. Similarly one can argue that in eq. (fT4|) 
only one type of B meson is involved and that it connects 
to only one of the 16 true Goldstone pions of the naive 
light quark theory (out of the total number of 256 pions). 
In other words both eq. (fT%)) and the $>b two-point cor- 
relator have the same normalization as in a theory with 
undoubled light and heavy quarks, such as continuum 
QCD. The one correlator where adjustment of normal- 
ization is required is the $„• — <!>„. correlator. Using naive 
fields as in ea. lfTS")) brings in an extra factor of 16 due to 



the trace over a 16 x 16 taste matrix and this factor of 
16 must be divided out. If one works with conventional 
staggered light quarks one would have an extra factor of 
4 rather than 16 relative to a pion correlator in a theory 
with undoubled fcrmions. 

To extract the matrix elements (tt| jjf\t)\B), the 3- 
point correlators must be fitted to: 

N w -1 N B -l 

cW(p„,p B ,t,T B )^ Yl E 

k=0 1=0 



(18) 



With this ansatz every second exponential (k or I odd) 
corresponds to an oscillatory (in time) contribution to the 
correlator, a characteristic feature of staggered fermions. 
We use Bayesian fitting methods (4f| and in most of our 
fits we kept N v = 1 and let Nb vary between 3 ~ 9. 
In order to avoid contamination from excited pions we 
dropped 5 to 8 points close to the pion source and made 
certain that within errors fit results did not depend on 
the number of points omitted. We have also tried some 
fits with N-x — 2 or 3. For some light quark masses and 
pion momenta one obtained results consistent with the 
N„ = 1 fits, in other cases, however, it was not possible 
to get very stable (with respect to Nb) fits. Hence, for 
our final results we rely on the N n — 1 fits. Since os- 
cillatory contributions are much more significant in the 
B channel than for pions, it was more important to al- 
low Nb to increase at fixed N n rather than the other 
way around. For priors in our Bayesian fits we use cen- 
tral values corresponding to energy splittings of about 
400MeV allowing, however for large ~100% widths. As 
priors for the amplitudes we allow for ranges typically 
between about —10 x A o and +10 x A o, where Aqo is 
the groundstate amplitude. Fit results for groundstate 
energies and amplitudes were very insensitive to choices 
for priors. In Figs. 5 & 6 we show examples of fit results 
for (jji(t)) for pion momentum -^2— (0,1,1) on one of 
the coarse ensembles. For comparison we show in Fig. 7 
results for a fit to the B correlator that was done simul- 
taneously with the fit to the 3-point correlator. 

The main goal in all our fits is to extract the ground 
state amplitudes A^ 00 which lead directly to the form 
factors f» and f± via eq.©. 



/II 



f± = 



A 00 (V ) 



A oa (V k ) 



(19) 



Here AqqIV^) includes contributions from all cur- 
rents with appropriate matching factors as dictated by 
eqns. (fl"2"| ) and (flU)) . ^ and <^b are the ground state am- 
plitudes from the pion and B meson correlators, respec- 
tively. Ctt is correctly normalized as in continuum QCD. 
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FIG. 5: Fit result for (J< 0) (t)) x e Mst x e E ' ( - T ' t) versus i for 
pion momentum (0, l,!)^ 2 -. The horizontal time axis has 
been rearranged so that the B meson source is at t/a = 1 and 
the pion source at t/a = T = 16. 
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FIG. 7: Fit results for the B meson correlator 



uoarrif 


uoam q 


p* = (000) 


(001) 


(011) 


(111) 


0.005 


0.005 


1.486(34) 


1.326(36) 


1.221(38) 


1.127(78) 


0.007 


0.007 


1.538(32) 


1.328(33) 


1.185(49) 


1.013(115) 


0.010 


0.005 


1.582(37) 


1.322(37) 


1.201(60) 


1.053(121) 


0.010 


0.010 


1.584(44) 


1.322(39) 


1.212(52) 


1.064(138) 


0.010 


0.020 


1.581(30) 


1.372(29) 


1.253(37) 


1.117(72) 


0.020 


0.020 


1.508(48) 


1.378(43) 


1.264(47) 


1.121(83) 



TABLE III: Results for the form factor /y in GeV 1/2 from 
coarse MILC lattices. Errors are statistical errors coming 
from Bayesian bootstrap fits. 
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V. CHIRAL EXTRAPOLATIONS 

The form factors listed in Tables III and IV are for an 
unphysical world with u and d quark masses larger than 
in reality. They need to be extrapolated to the physical 
chiral limit. In ref.[17| linear extrapolations were per- 
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FIG. 6: Same as Fig.5 but for (J< 0) (t)} 



Results for the form factors /m and f± for the different 
ensembles and pion momenta are summarized in Tables 
III & IV. 



uoarrif 


uoam q 


P* = (001) 


(on) 


(111) 


0.005 


0.005 


1.543(205) 


0.703(44) 


0.515(39) 


0.007 


0.007 


1.082(31) 


0.606(31) 


0.433(48) 


0.010 


0.005 


1.128(37) 


0.662(36) 


0.413(38) 


0.010 


0.010 


1.235(90) 


0.657(37) 


0.449(43) 


0.010 


0.020 


1.029(21) 


0.611(19) 


0.423(23) 


0.020 


0.020 


1.097(140) 


0.597(24) 


0.397(20) 



TABLE IV: 
GeV" 1/2 . 



Same as Table III for the form factor f± in 
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formed on the three points with uoarrif = 0.01, the only 
simulation data available at that time. Here we carry out 
this important step in several different ways : 
1) linear extrapolation with only the full QCD (nif = 
m q ) data, 2) staggered chiral perturbation theory 
(SChPT) with full QCD data, 3) continuum ChPT with 
full QCD data, and 4) SChPT simultaneously to both 
full QCD and partially quenched data. We use the last, 
most involved chiral extrapolation for our final answer, 
but make certain that the other methods give agreement 
within quoted errors. 

The formulas of chiral perturbation theory for heavy- 
light form factors have the general form, 

f\\/± = c [l + Sf\i/ ± + cim q + c 2 (2m f + m s ) + ....]. (20) 

<5/|| and dfx are the chiral log terms and they have 
been determined at lowest order in 1/M for continuum, 
quenched and partially quenched QC D liq. 14711 and also 
specifically for staggered light quarks [34[ [33, [36| . We use 
the formulas of Aubin & Bernard [36[ that include "taste 
symmetry breaking" lattice artifact terms that come in 
at (D(a 2 ). Sfu and Sf± are functions of E n (more gener- 
ally of v -p n , with equal to the B four velocity and 
the pion momentum) and also depend on the B* Bit cou- 
pling gBn- A possible source of concern in using ChPT 
may be the presence of pions with E v > 2m n in our sys- 
tem. Higher order terms in continuum heavy-light ChPT 
are discussed, for instance, in ref . [48j| . There it is argued 
that for some processes, effects of higher order terms can 
be taken into account by allowing for a correction to gs^ 
linear in E v . In our fits to SChPT and to continuum 
ChPT formulas we will let gs-w vary as one goes from 
one E„ value to another, together with cq, c\ and C2. 

The discussion in the previous paragraph shows that 
chiral extrapolations are most conveniently carried out at 
fixed values of E v . In order to do so, one needs to inter- 
polate the data of Tables III & IV to fixed common values 
of E-x for each light quark mass. We have explored several 
different ways to perform the interpolations. We take ad- 
vantage of various ansaetze that have been developed in 
the literature to model the q 2 dependence of form factors 
f+(q 2 ) and fo(q ). A ppe ndix A summarizes commonly 
used ansaetze [3 HJ M [M IM ED- IM • The formulas 
given there refer to form factors in the real world, i.e. 
when using physical values for pion, B and B* masses. 
We use them here also as guides for sensible interpo- 
lations even at unphysical meson masses. We caution, 
however, that quantities such as q^ ax = (Mb — m-*) 2 , 
M B , — q^ax 01 the B ir threshold (Mb + va^) 2 are sensi- 
tive to the pion mass. A particular ansatz may not work 
as well with unphysical masses compared to in the chiral 
limit, or vice versa. For instance, we find that away from 
the chiral limit, i.e. when we are interpolating the data 
of Tables III and IV, the most popular three parameter 
Becirevic-Kaidalov (BK) ansatz [49( often fails to accom- 
modate our data. After the chiral extrapolation the BK 
ansatz works well, however. On the other hand we were 
often (though not always) able to fit data away from the 




q 2 (GeV 2 ) 

FIG. 8: Comparison between a 3 parameter BK fit and a 4 
parameter BZ fit to uoam q = uoarrif = 0.01 data. 



chiral limit using just a simple single pole ansatz for both 
/o and /+, but as expected, a single pole ansatz does not 
work for f+(q 2 ) in the physical limit. The ansatz that 
works best both away and at the chiral limit is one due 
to Ball & Zwicky (BZ)[5(|, which actually is the same as 
the four parameter BK parametrization (see Appendix). 
Examples of fits used in the interpolations are given in 
Figs. 8 & 9. So, for each (m/,m g ) we first convert the 
results for /|| and f± of Tables III & IV to fo and /+, 
interpolate using the BZ ansatz and finally convert back 
to /[| and fx at fixed values of E„. These are then ex- 
trapolated to the chiral limit using the SChPT formulas 
described above. One advantage of carrying out interpo- 
lations for /o and /+ rather than directly in fu and fx, 
is that the kinematic constraint /o(0) = /+(0) is easily 
incorporated into the ansaetze and one can use the fo 
data to constrain the normalization of /+ as well. 

In Figs. 10 & 11 we show chiral extrapolations for /y 
and fx respectively for several values of E„ . We compare 
results from simple linear extrapolations of full QCD data 
and from using the SChPT formulas simultaneously for 
both full QCD and partially quenched QCD data. One 
sees that the difference between the two ways of perform- 
ing the chiral extrapolations is small. For f» the central 
values at the chiral limit differ only by 0.1 ~ 2.4%. For 
fx the differences are again < 2.4% for E^ > 0.7GeV. 
For E« < OJGeV (not shown in Fig.9) linear and SChPT 
extrapolations of fx can differ by 4 ~ 6%. On the other 
hand this is the kinematic region where statistical (and 
interpolation) errors are large for some combinations of 
m f and m q (see third column in Table IV) and it is not 
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FIG. 9: Same as Fig. 8 for uoam q = 0.02, uoarrif = 0.01. 



possible to disentangle statistical and chiral extrapola- 
tion errors. f± at small E v translates into /+ at large 
q 2 . One can speculate that the large statistical errors we 
are finding in this kinematic region may be due to the 
proximity of the Mb* pole. This, however, for reasons 
we do not understand, was not evident in the partially 
quenched data of reference [ItJ ■ 

For the SChPT extrapolations, as mentioned above, 
we have let qb-k float and be one of the fit parameters. 
We find that the "effective" g Bv ranges between 0.0 and 

0. 2 and decreases with E n , although within large errors. 
In Fig. 12 we show a comparison of fits with and with- 
out the 0{a 2 ) corrections in the ChPT theory formulas, 

1. e. a comparison between continuum and staggered full 
QCD ChPT. Again the differences are very small. In 
summary, for most of our data points, different ways of 
carrying out the chiral extrapolation, including using no 
input from ChPT at all, lead to a spread in extrapolated 
results of only 2.5% or less. This indicates that contri- 
butions we have neglected in the ChPT formulas, such 
as 1/M corrections, higher order (in pion momentum) 
terms, finite volume effects etc. are not important. We 
take as central values for the form factors /|| & f± in the 
chiral limit, the results coming from the SChPT extrapo- 
lation using both full QCD and partially quenched data. 
The combined statistical and chiral extrapolation errors 
are discussed in the next section. 

In Fig. 13 we compare results from one of the MILC fine 
lattice ensembles with the coarse lattice data discussed 
so far. For this comparison both the coarse and fine 
data points include only the (J^) contributions through 
0(a s ), i.e. the first lines in eqns. (fT2") ) or (fTB")) respectively. 
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FIG. 10: Chiral extrapolation of /y . The circles to the left 
show the extrapolated values at the physical limit. 



1.5 



> 

o 





0.5 



— Staggered Chiral perturbation theory fit 

— linear fit 

□ E s =0.8GeV 

E^l.O 
A E n =l.l 
< E =1.2 




0.2 



0.4 
m Im 

q s 



0.6 



0.8 



FIG. 11: Same as Fig.8 for f ± . 



Since we have shown that the higher order currents have 
minimal effect, we believe meaningful scaling tests can be 
carried out with just (jji }. For which is the main 
contributor to the phenomenologically relevant form fac- 
tor f + (q 2 ) , one sees that the fine lattice point falls nicely 
on the fixed E w curve determined by coarse lattice data. 
Our statistical error is large currently for the /|| point 
from the fine lattice. Within these large errors there is 
consistency between the coarse and fine lattices for /|| 
as well. We conclude from this exercise that there are 
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FIG. 12: Chiral extrapolations with (staggered ChPT) and 
without (continuum ChPT) 0(a 2 ) corrections. 



1.5 



0.5 



6 



f^GeV"" 2 ) 
O fine lattice 
O UGeV 1 ' 2 ) 



□ 



0.1 



0.2 



0.3 
m /m 

q s 



0.4 



0.5 



FIG. 13: Comparison of coarse lattice data with some results 
from one of the fine MILC ensembles. Shown are results for 
fx and /n at E n = 0.79GeV. 



no indications of large discretization effects in the form 
factor calculations on MILC coarse lattices. Such errors 
are smaller than current statistical errors. Eventually it 
would be desirable to carry out a more thorough scaling 
test, once more data on fine lattices at several light sea 
quark masses become available. 



VI. RESULTS FOR FORM FACTORS / 4 
AND f (q 2 ) IN THE CHIRAL LIMIT 



and 
in the 



We convert the chirally extrapolated fn(E n ) 
f±(E 7T ) to the form factors f+{q 2 ) and /o(<? 2 ) 
physical limit. These are shown in Fig. 14 and tabu- 
lated in Table V. For comparison we also plot in Fig. 14 
the data presented in ref. [17(. One sees that changes 
are minimal in spite of all the improvements included in 
our new results. This indicates that the approximations 
that were made previously and that we are systemati- 
cally improving upon, such as partial quenching, linear 
chiral extrapolations, working with currents at lowest or- 
der in 1/M, did not drastically affect the theory. The 
solid curves in Fig. 14 are fits to our new results using the 
Ball-Zwicky (BZ) [5(| parametrization of / + and / . We 
have also tried fits to other parametrizations, described 
in the Appendix, including the Becircvic-Kaidalov (BK) 
49 1, Richard Hill (RH) [20| and a series expansion (SE) 
18 , [U [5l|, parametrization. The RH parametriza- 
tion fit is essentially indistinguishable from the BZ fit. 
The BK fit is also a good fit to our data although not 
quite as good as the first two. This should not be surpris- 
ing, since the BK fit has only three parameters to tune 
whereas the BZ and RH fits are both four parameter 
fits. Any further parameters, however, are very poorly 
determined and do not help in the fit. Another class of 
fit ansaetze, the series expansion (SE) fits, are discussed 
in the Appendix. The main reason we are interested in 
obtaining a good analytic parametrization of the form 
factors, is to facilitate partial integration of differential 
decay rates, as discussed below. These parametrizations 
can also be used to try and extrapolate to lower q 2 where 
lattice data are currently not available. 

The statistical plus chiral extrapolation errors for 
f+(q 2 ) lie between 7 ~ 10% depending on q 2 . They are 
smaller for the form factor fo{q 2 )- For q 2 > 16GeV 2 , the 
range we will be focusing on, the average error for f + (q 2 ) 
comes out to be ~ 8%. In Table VI we list this average 
statistical plus chiral extrapolation error together with 
estimates of systematic errors from other sources. These 
other systematic errors are dominated by the ~ 9% un- 
certainty in higher order matching of the heavy- light cur- 
rents. 

The differential partial decay rate for B — > tt Iv, ignor- 
ing the charged lepton mass, is given by, 



dT 
dq^ 



^pIWuA 2 \f + (i 2 )\ 2 



(21) 



where Gf is the Fermi constant and p v the magnitude 
of the pion three momentum in the B rest frame. Know- 
ing f+{q 2 ) then allows us to evaluate ry~p ^2 and also 
integrate this quantity over different q 2 bins. We take 
our best fit, the BZ fit shown in Fig. 14, and integrate to 
obtain, 
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FIG. 14: Form factors f+(q 2 ) and foil 2 ) in the chiral limit. 
The black squares and triangles are the new and final results 
for f± and /o respectively. For comparison, the data from 
ref.[17l] are also shown as circles. The full black curves follow 
a BZ parametrization fit (see text) to the new data. Errors 
are combined statistical and chiral extrapolation errors. 



source of error 


size of error (%) 


statistics + chiral extrapolations 


8 


two-loop matching 


9 


discretization 


3 


relativistic 


1 


Total 


12 



TABLE VI: Estimate of percentage errors in f+(q 2 ) for q 2 > 
16GeV 2 . 



Fit 


f+(q 2 = o) 


/ dg 5 '/ 


V ub \ 2 [ps- 1 ] 






< q 2 < q'Lax 


16GeV 2 < q 2 < q 2 max 


BZ 


0.27(2) (4) 


6.00(96) (1.68) 


1.46(23)(27) 


BK 


0.26(2)(4) 


6.03(96) (1.69) 


1.31(21)(25) 


RH 


0.27(2) (4) 


5.99(96) (1.68) 


1.45(23) (27) 



TABLE VII: Partially integrated differential decay rates and 
/+(0) using several parametrizations. The first error reflects 
statistical and chiral extrapolation uncertainties. The second 
error is due to remaining systematic errors. 



VII. ESTIMATING \V U 



q 2 [GeV 2 ] 


U(i 2 ) 


Mi 2 ) 


15.23 


0.649(63) 


0.475(26) 


16.28 


0.727(64) 


0.508(25) 


17.34 


0.815(65) 


0.527(25) 


18.39 


0.944(66) 


0.568(24) 


19.45 


1.098(67) 


0.610(24) 


20.51 


1.248(97) 


0.651(25) 


21.56 


1.554(156) 


0.703(26) 



TABLE V: Form factors f+(q 2 ) and fo(q 2 ) in the chiral limit. 
Errors shown are combined statistical and chiral extrapolation 
errors. 



The first error is the combined statistical plus chiral ex- 
trapolation error and the second the sum of all other sys- 
tematic errors added in quadrature. Eq. (122p is the main 
result of this article. It serves as basis for determinations 
of the CKM matrix element Similar integrals us- 

ing other parametrizations and over other q 2 ranges are 
summarized in Table VII together with /+(0) from the 
different fits. The second error for /+(0) and for the inte- 
grated rate over the entire q 2 range includes an additional 
10% systematic uncertainty in f + (q 2 ) which is not part 
of Table VI and which comes from the extrapolation into 
the low q 2 region. 



In this section we combine the lattice results of the 
previous section with experimental input for B — > tt Iv 
branching fractions and extract an estimate for \V u b\- In 
order to do so, we use results from the Heavy Flavor 
Averaging Group (HFAG) [53| and rely on its analysis 
of the current experimental uncertainties. The HFAG 
gives preliminary averages of BaBar, Belle and CLEO 
results (as of the conferences of Summer 2005) for the 
integrated branching fraction B{B° — > 7r~ They 
quote [1.35 ± 0.08 ± 0.08] x 10~ 4 for < q 2 < q 2 max and 
[0.40 ±0.04 ±0.04] x 10~ 4 for q 2 > 16GeV 2 . Combining 
this with eq. (|22|) and a B° lifetime of 1.536 ps [54| leads 
to, 

[K&l = 4.22(30)(51) x 10~ 3 , q 2 > IQGeV 2 (23) 

where the first error is experimental (7%) and the sec- 
ond is the total lattice error (12%). The result (|23| is 
consistent, at the one a level, with the preliminary value 
3.78(25)(52) x 10" 3 obtained by the Fermilab/MILC col- 
laboration using the same HFAG branching fraction av- 
erages [55| . 

VIII. SUMMARY 

We have completed a determination of the B meson 
semileptonic form factors f+(q 2 ) and fo{q 2 ) using state- 
of-the-art Lattice QCD methods. Our calculations em- 
ploy unquenched gauge configurations, created by the 
MILC collaboration, that incorporate vacuum polariza- 
tion effects from two very light flavors plus strange sea 
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quarks. Both the sea and the valence light quarks are 
simulated using a highly improved staggered quark ac- 
tion. This action allows us to work close enough to the 
chiral limit, so that chiral extrapolations to physical pi- 
ons are mild and do not introduce large uncertainties. 
Our results for f+(q 2 ) can be combined with experimen- 
tal branching fraction data to extract the CKM matrix 
element \V u b\- This quantity is a crucial ingredient in 
tests of the unitarity triangle and in solidifying our un- 
derstanding of CP violation in the Standard Model. 

The total lattice error in the f+(q 2 ) form factor pre- 
sented here, and hence also in |14f,|, is at the ~ 12% level. 
This error is dominated by uncertainty in higher order 
perturbative matching of heavy- light currents and by sta- 
tistical errors. One of the goals of the HPQCD collabora- 
tion is to carry out higher order matching for heavy- light, 
light-light and heavy-heavy currents and four-fermion op- 
erators 56] . It has become increasingly obvious that such 
calculations are necessary for accurate lattice determina- 
tions of form factors, decay constants and mixing pa- 
rameters at the ~ 5% level. For form factor calculations, 
more work is also required to reduce statistical errors. 
One approach that helped significantly in the HPQCD 
collaboration Jb determination [37} is to explore differ- 
ent smearings of sources and sinks in correlators. Wc 
plan to investigate this in the future. 

Lattice results for f+(q 2 ) exist at the moment only for 
q 2 > lQGeV 2 and this forces us to use only part of the 
available experimental branching fraction data. In or- 
der to take advantage of all the experimental data and 
thereby reduce experimental errors in | V u b | , other meth- 
ods such as sum rules approaches are currently employed 
[50L |j57j to cover the low q 2 region. A lattice approach 
to f+(q 2 ) determinations at low q 2 , "Moving NRQCD", 
was introduced many years ago and work is in progress 
to implement this method using our highly improved 
quark and gauge actions [H, [59[ . We can look forward 
to the next major improvement in lattice determinations 
of \V u b\ once two- loop matching, reduction in statisti- 
cal errors and control over the entire q 2 range has been 
achieved. 



APPENDIX A: PARAMETRIZATION OF FORM 
FACTORS 

Most parametrizations start from a dispersive repre- 
sentation of the form factors. 



Mi 2 ) = 



fo(q 2 



1 



(1-9 2 ) 



dt 



*+ 

r-2 



Im[f+(t)] 
t — q 2 — ie 



(1-?) (1-af) 
r\ + r 2 



(1 - q*/0) 



(Al) 



(A2) 



where q 2 = q 2 /M Bte and t + = (M B + m^) 2 . The kine- 
matic constraint /+(0) = /o(0) is automatically satis- 
fied. This is a 4 parameter parametrization of /+ and /o, 
sometimes called the "4 Parameter BK Parametrization" 
(Becirevic & Kaidalov) j4§] . Examples of generalizations 
and special cases developed in the literature are given 
below. They differ only in the parametrization of /+. 



3 Parameter Becirevic-Kaidalov (BK) [49( 
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fo(q 2 ) = 
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4 Parameter Ball-Zwicky (BZ) [5C 
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^ _ /+(0) 
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4 Parameter R.Hill (RH) |2C 



/+(«f) = 



2 . f+(0)(l-S-q 2 ) 
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Another class of parametrizations, advocated in (l9l . [BTI . 
[52} and also discussed in [2l[ , is based on a series expan- 
sion of f+{q 2 ) around some q 2 = t . For better conver- 
gence of the series it is customary to convert to a new 
variable z(g 2 ,io), where following [n| [2l| we take, 

z{q »*o) = , = — - ■ (A7) 
y/t+ -q 2 + - t Q 

We refer to the literature for further discussion of the 
merits of this transformation. The form factor f + (q 2 ) 
can then be expanded as a power series in z(q 2 ,to), 

k 
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to / Qmax 


J llq 1 1 


V ub \ 2 [ps- 1 ] 




< q 2 < q^nax 


16GeV 2 < q 2 < q 2 max 


0.0 


5.78 


1.32 


0.2 


5.78 


1.33 


0.4 


5.79 


1.33 


0.6 


5.80 


1.33 


0.8 


5.79 


1.33 


1.0 


5.79 


1.33 



TABLE VIII: Partially integrated differential decay rates us- 
ing the SE parametrization for different choices of to- No error 
estimates have been made. 



The "Blaschke" factor P(q 2 ) must take into account any- 
isolated poles below the Bn threshold at q 2 = t + . Wc 



set P(q 2 ) = z(q 2 ,Mg,) to take care of the B* pole. For 
<&(<7 2 , we take the expression given in ref.[19| (with 
simplified X )■ We combine the ansatz (A8) for f+{q 2 ) 
with k m ax = 2 together with (A4) for fo{q 2 ) to get 
another (and our last) 4 parameter ansatz, 



Series Expansion (SE) 



a + aiz{q 2 , t a ) + a 2 z 2 
z(g 2 ,Ml,)<% 2 ,*o) 



(A9) 



We have explored the SE parametrization but not be- 
cause we needed a better analytic expression to cover 
the range q 2 > 16GeV 2 of our lattice results. Rather, 
we did so to assess how reliable is the information on 
the shape of the form factors at lower q 2 that we are 
getting from the BK/BZ/RH type parametrizations. We 
find that good fits to the simulation data can be obtained 
with the SE parametrizations as well and that they are 
insensitive to the value of to. Table VIII gives results 
for integrated differential decay rates. Not surprisingly, 
agreement is found with Table VII for q 2 > 16GeV 2 . Re- 
sults are systematically slightly lower than in Table VII 
if one includes the entire q 2 range. One can summarize 
by saying that we do not have evidence for any strong 
dependence on the choice of ansatz employed to extrap- 
olate lattice data to lower q 2 values. Nevertheless, we 
believe it is important to look for additional information 
on the form factors at q 2 < 16GeV 2 , either from experi- 
ment or from theoretical models, if one wants to discuss 
the entire q 2 range at the present time. This has already 
been done in the recent literature [l|| [2l|, [5(| • Alterna- 
tively, getting lattice results directly in the low q 2 regime 
using for instance Moving NRQCD would also solve this 
problem. 



APPENDIX B: ERRATUM: B MESON 
SEMILEPTONIC FORM FACTORS FROM 
UNQUENCHED LATTICE QCD [PHYS.REV.D 
73, 074502 (2006)] 

Due to a normalization error in one of our analysis codes 
many results for the form factor f± were evaluated in- 
correctly. In particular, Table IV should be replaced by, 







f± [Gev 


1/2] 




uoanif 


uoam q 


V, = (001) 


(on) 


(HI) 


0.005 


0.005 


1.543(205) 


0.994(62) 


0.892(68) 


0.007 


0.007 


1.082(31) 


0.857(44) 


0.750(83) 


0.010 


0.005 


1.128(37) 


0.936(51) 


0.715(66) 


0.010 


0.010 


1.235(90) 


0.929(52) 


0.778(74) 


0.010 


0.020 


1.029(21) 


0.864(27) 


0.733(40) 


0.020 


0.020 


1.097(140) 


0.844(34) 


0.688(35) 



The corrected form factors f+(q 2 ) and fo(q 2 ) in the chiral 
limit become (previous Table V): 



q 2 [GeV 2 ] 


/+(<z 2 ) 


fo(q 2 ) 


17.34 


1.101(53) 


0.561(26) 


18.39 


1.273(99) 


0.600(21) 


19.45 


1.458(142) 


0.639(23) 


20.51 


1.627(185) 


0.676(41) 


21.56 


1.816(126) 


0.714(56) 



We restrict the range in q 2 to values where simulation 
data before chiral extrapolation are available for all light 
quark masses and omit two lower q 2 points that were 
obtained previously through small extrapolations in E v . 

The new error budget is given by (previous Table VI): 



source of error 


size of error (%) 


statistics + chiral extrapolations 


10 


two-loop matching 


9 


discretization 


3 


relativistic 


1 


Total 


14 



and the partially integrated differential decay rates be- 
come (previous Table VII): 



Fit 


U(q 2 =0) 


/ Iq 1 / 


V ub \ 2 [ps" 1 ] 






< q 2 < q 2 nax 


16GeV 2 <q 2 < q 2 max 


BZ 
BK 
SE 


0.31(5)(4) 
0.31(5)(4) 
0.30(5)(4) 


9.10(1. 82)(2. 55) 
9.30(1.86) (2.60) 
9.35(1.87) (2.62) 


2.07(41)(39) 
2.13(43)(40) 
2.02(40)(38) 



Equation (22) should be replaced by, 
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/ -73 ^=2.07(41X39) jw" 1 , 

which signifies an increase on our previous incorrect value 
( 1.46(23)(27) ps^ 1 ) by 1.7 times the previously quoted 
total error. 

The new eq.(22) leads to a new eq.(23), 

\V ub \ = 3.55(25)(50) x 10~ 3 , q 2 > 16GeV 2 . 



Accordingly, the last two sentences in the Abstract 
should be changed to, 



We calculate the form factors f+(q 2 ) and fo(q 2 ) in the 
chiral limit for the range 16 GeV 2 < q 2 < q^aj. and ob- 
tain jf^ y2 [dT/dq 2 ]dq 2 I \V uh \ 2 = 2.07(57) ps-\ Com- 
bining this with a preliminary average by the Heavy 
Flavor Averaging Group (HFAG'05) of recent branch- 
ing fraction data for exclusive B semileptonic decays 
from the BaBar, Belle and CLEO collaborations, leads 
to |K&| = 3.55(25)(50) x 1(T 3 ." 
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